Turbulent fronts of quantum detonation in molecular magnets 
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Dipolar-controlled quantum deflagration going over into quantum detonation in the elongated 
Mni2 Ac molecular magnet in a strong transverse field has been considered within the full model. 
It is shown that within the dipolar window around tunneling resonances the deflagration front is 
non-flat. With increasing bias, dipolar instability makes the front turbulent, while its speed reaches 
sonic values, that is a signature of detonation. 
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Magnetic burning or deflagration, similar to the usual 
chemical burning, has been experimentally observed on 
elongated crystals of molecular magnet (MM) Mni2 Ac in 
Refs. [HE]. Deflagration is controlled by the temperature 
rising as the result of the relaxation from the metastable 
to stable state of this bistable system (burning), accom- 
panied by the heat conduction toward the cold unburned 
region. This leads to the formation of a flat burning front 
moving at a constant speed controlled by the energy re- 
lease set by the bias magnetic field Bz applied along the 
easy axis of magnetic molecules. This discovery triggered 
theoretical [3j and further experimental [4, 5j research. In 
these experiments no transverse field was applied and the 
front speed was as low as 1-15 m/s. Fast moving burning 
fronts in Mni2 Ac initiated by a fast magnetic field sweep 
have been observed in Ref. [6j. In this region, defiagra- 
tion can go over into detonation with a speed comparable 
to the speed of sound [7J. Crossover from relaxation to 
defiagration in Mni2 Ac in a strong transverse field has 
been experimentally studied in Ref. [8j. 

Since MM are famous exponents of the resonant spin 
tunneling [9l |10] (see Ref. pTj for a review) , manifesta- 
tions of the latter were expected in magnetic defiagra- 
tion. Indeed, maxima of the front speed at tunneling 
resonances were seen in Refs. [11 [4]. The simplest way to 
explain these maxima is to use, within the standard de- 
fiagration theory, the Arrhenius relaxation rate V{Bz-,T) 
having maxima at Bz corresponding to tunneling reso- 
nances [U [3] . 

It was further suggested that in the case of a strong 
tunneling directly out of the metastable ground state, 
caused by a strong transverse magnetic field, fronts of 
non-thermal quantum or "cold" defiagration are possi- 
ble [121 US]- Here, instead of the temperature, tunnel- 
ing is controlled by the dipolar field of the crystal that 
can block and unblock tunneling by setting magnetic 
molecules on or off resonance. Within the simplified Id 
approximation it was shown that the magnetization in 
the moving front is adjusting self-consistently, so that 
the dipolar field unblocks tunneling within the front core. 
The theory of quantum defiagration has been generalized 
[T4[ ITS] to include both quantum and thermal effects via 
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Figure 1: (Color online). Dipolar instability of a smooth front 
of tunneling. The leading parts of the front (center) are helped 
by the lagging neighbors while hampering the motion of the 
latter. 

the relaxation rate T{Bz-,T) numerically calculated from 
the density-matrix equation [16j. The main finding was a 
strong quantum acceleration of the front within the dipo- 
lar window around the resonance values of the external 
field Bz = Bk: 

Bk-Bf-^ <Bz<Bk^Bi''\ (1) 

where B^^^ =52.6 mT is the dipolar field produced by 
the fully magnetized long crystal and B^J^^^ = 72.9 mT 
[15j.The front speed was supersonic and apparently di- 
verging towards the right end of the dipolar window, see 
Fig. 15 of Ref. fT5^ and discussion therein. Indeed, as the 
dipole-dipole interaction (DDI) is instantaneous, there is 
no limitation on the front speed. 

While fast quantum defiagration in MM is avaiting ex- 
perimental confirmation, a theoretical problem remains. 
The Id approximation assumes that the magnetization 
and dipolar field in the front depend only on the coordi- 
nate z along the sample (that coincides with the MM easy 
axis). This is a crude approximation, and there should 
be a dependence on x and y that renders the front non- 
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flat. Moreover, dipolar- controlled front should be unsta- 
ble at a small scale because portions of the front that are 
ahead of the neighbors create dipolar field on the latter 
that hampers their motion, while the lagging neighbors 
create a field that accelerates the leader (see Fig. [T]). 
This instability refiects the fact that a long MM crystal 
tends to split up into domains parallel rather than per- 
pendicular to the z axis [17]. As the result, the front can 
become turbulent with a chaotic dipolar bias rendering 
many magnetic molecules inside the front core off reso- 
nance. The question is whether the front speed drops to 
that of the regular slow burning jlj |18] or remains fast. 
In the latter case one has to speak of quantum detonation 
rather than of quantum defiagration. 

The purpose of this Letter is to present a full 3d theory 
of fronts of tunneling in molecular magnets and clarify 
the role of the dipolar instability. A 3d model with DDI 
presents a numerical challenge, still its solution is possible 
on modern workstations. 

The system of equation describing defiagration with 
quantum effects in molecular magnets [14j includes the 
rate equation describing relaxation of the metastable 
population n (1 < n < 1) 



n(t,r) = -r[5tot,.(r),T(r)] n(t,r) 



,(eq) 



(2) 



Here the relaxation rate F depends on the total bias field 

5tot,.(r) = + = B,^{Sg^B/vo)D,,{r), (3) 

(external plus dipolar field) and the temperature at a 
given point, while n^^^\T) is the equilibrium metastable 
population (set to zero below). Here vq = o?c is the unit- 
cell volume, a and c are lattice spacings, and Dzzij^) is 
the reduced dipolar field. Another equation is the heat 
conduction equation that can be conveniently written in 
terms of the thermal energy £ per magnetic molecule 



8(t, r) = V • nV£{t, r) - /\Eh{t, r). 



(4) 



Here = 2SgiiBBz in the source term of this equation 
is the released energy per molecule, k is thermal diffusiv- 
ity (a crude estimate n 10~^ m^/s [HIS]). £ and T are 
related via the measured [19j heat capacity C = d£/dT. 

While the Eqs. |2]and|4] are the same as the standard 
defiagration equations [3j, the resonance form of F and 
its crucial dependence on the dipolar field B^^"^ (r) that is 
defined by n(t, r) makes a big difference. For the realistic 
model of Mni2 Ac with the uniaxial part of the Hamilto- 
nian —DS^ — AS^^ tunneling resonances are achieved at 
Btot,z Bkm, where 



gi-iBBkm = A: + (m^ + (m + kf) A] 



(5) 



k = 0,1,2,... and m = —S^—S + 1,... label 
the metastable spin states. The small fourth-order 
anisotropy A splits k resonances into m multiplets. With 
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Figure 2: (Color online). Relaxation rate of Mni2Ac vs 
temperature and longitudinal magnetic field in the trans- 
verse field B± = 3.5 T. One can see the ground-state res- 
onance at Bz = 0.522T and the first-excited- state resonance 
at Bz = 0.490 T for /c = 1 multiplet. 



Mni2 Ac parameters from Ref. [20j and using the density- 
matrix-equation method of Ref. [16j, for the transverse 
field B± = 3.5 T one obtains F around the first reso- 
nance, /c = 1, shown in Fig. [2j At such strong B± only 
the resonances with m = — 10 (ground-state resonance) 
and m = — 9 are seen, while all the other broadened away. 

The most challenging part of the work is numerical 
calculation of the dipolar field produced by the sample. 
Dzzij^) in Eq. (|3| is defined by the lattice sum Dzzij^i) = 



^ . (j){Yj - Yi)Gz{Yj), where 
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To calculate this sum, one can introduce a small sphere 

1 /3 

of radius ro around satisfying Vq <C tq <C where 
L is the (macrocopic) size of the sample. Assuming that 
^■^(r) does not change at the scale ro, one obtains 



(r) = - / dv'^(v' - v)cjz(v') + cjz(v)D^^t\ 

^0 J|r'-r|>ro 



(7) 

where v = 2\s the number of magnetic molecules per unit 
cell ^0 for the body-centered tetragonal Mni2 Ac and the 
lattice-dependent D^zz^^ = 2.155 comes from summation 
over the small sphere. The restriction on the integration 
can be removed by subtracting o-z{r) from the integrand 
and adding the counter-term. Finally one obtains 

Dzz{r) = - [ dv'c^iv' -v){cjz(v')-cjz(v)) 
vo J 

+ az{r){uVzz{r)-kD), (8) 



where ko = ^"kv /3 — D 
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14.6. Vzz{^) is the reduced 



magnetostatic field given by an integral over the sample 
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5, = 0.56 T 
v= 1674 m/s 



Figure 3: (Color online). Front speed within the 3d model for 
a strong transverse field [B± — 3.5 T) in the vicinity of the 
ground-state tunneling resonance at Bz — — 0.522 T. 



Figure 4: (Color online). Profile of the metastable population 
n in the ?>d model of fronts of tunneling in a 0.2 mm Mni2Ac 
crystal at B^ = 3.5 T and B^ = 0.5 T (upper) and 0.56 T 
(lower). 



surface. For a rectangular sample with dimensions 2Lx x 
2Ly X 2Lz and —Lx < x < Lx etc. one has [21j 

(r) = arctan 

Vx,Vy,Vz=±l 

{Lx + Vxxy^ {Ly + r]yy) (L^ + iq^z) 
\J {Lx + r]xxf + {Ly + r]yyf + (L^ + r]zzf 

-\-{x^y)^ in total 16 different arctan terms. Inside 
a uniformly magnetized long Mni2 Ac crystal Eq. (|8| 
yields Dzz = 10.53 that in real units corresponds to 
B^^^ = 52.6 mT. The same value of the dipolar field 
follows from the experiment |22| . 

For numerical solution, box-shape samples with dimen- 
sions Lx X Ly X Lz and ^ < x < Lx etc. are considered, 
and the problem is discretized with a rectangular grid. 
This makes the integral in Eq. (|8| a sum that can be cal- 
culated using the fast Fourier transform (FFT). In this 
work, numerical solution was implemented in Wolfram 
Mathematica, whereas the dipolar sum was computed 
using FFT based and very fast List Convolve command. 
After discretization of the Laplace operator in Eq. Q the 
system of equations for quantum deflagration becomes a 
system of nonlinear ordinary differential equations. It has 
been solved by the fixed-step compiled 4th-order Runge- 
Kutta method and, preferred, by somewhat more effi- 
cient 5th-order Butcher's Runge-Kutta method making 6 
evaluations per step. Unfortunately, long-range DDI pre- 
vents parallelization of the code. Instead oi 8^ it is more 
convenient to use £ / + n as one of unknown func- 
tions. The grid sizes were 64 x 64 x 300, 100 x 100 x 200, 
and 100 x 100 x 300. To reduce the amount of calcula- 
tions, symmetry conditions within xy planes have been 
imposed. 

Deflagration was ignited in the uniformly magnetized 
state n = 1 at the left end of the crystal, z = 0, by a 



quick temperature rise at the boundary. At higher 5^, 
deeper into the dipolar window, the resonance condition 
is fulfilled at a surface close to the left end, so that quan- 
tum deflagration begins spontaneously there. The right 
end of the sample was protected from spontaneous igni- 
tion by adding a semi-infinite fictituous sample with a 
fixed "down" magnetization, n = 1. The sample was con- 
sidered as thermally insulated, the initial temperature 
being almost zero. The front speed was determined from 
the time needed for the front to arrive at the right end 
or from the time dependence of the average metastable 
population. Computations were done for Lx = Ly = 0.2 
mm [4J and 1 mm [2\ . One burning event required about 
one day computations on a workstation. 

The results for the front speed in the vicinity of the first 
tunneling resonance {k = 1) are shown in Fig. [sj Beyond 
the dipolar window, Bz < 0.45 T and Bz > 0.573 T, there 
is no tunneling and the usual temperature-driven defla- 
gration with a smooth flat front and v ~ 50 m/s takes 
place. Entering the dipolar window from the left leads to 
a gradual increase of the effect of tunneling. The front 
speed increases and reaches sonic values, then the front 
becomes so fast that numerical calculations require a too 
long sample and become too difficult. Theoretically, the 
front speed diverges at the right end of the dipolar win- 
dow Bk + Bi^^ = 0.522 + 0.0526 = 0.573 T, where Bj, 
is the ground-state tunneling resonance and /c = 1. Al- 
though the front speed in the 3d model is lower than 
in the Id model, qualitatively the results are similar (cf. 
Fig. 15 of Ref. [15j). Here one cannot see any contri- 
bution of the excited-state tunneling resonance (peak at 
Bz = 0.490 T in Fig. |2|, the whole tunneling effect being 
due the the ground-state tunneling. 

The structure of the front of tunneling is shown in Fig. 
|4] Here n is represented as a function of x and z at 
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Figure 5: (Color online). Profile of the metastable population 
n in the 3d model of quantum deflagration in a 1 mm Mni2Ac 
crystal at B± = 3.5 T and = 0.55 T. 

y = 0. Color coding is blue for the cold unburned re- 
gions (n ~ 1) and red for hot burned regions (n ~ 0). 
Dipolar effect makes the front progressively non-flat, its 
central part leading as the speed is increasing, as seen in 
the upper image at 5^ = 0.50 T. Also with increasing Bz 
the front becomes gradually non-smooth. However, ther- 
mal burning tends to smoothen the roughness created 
by the dipolar instability discussed above as large tem- 
perature gradients lead to increased heat conduction and 
equilibration along the front. Whereas in the absence of 
thermal effect (cold deflagration) there is a spatially ir- 
regular unburned metastable population behind the front 
[T2l [T3]. here a complete burning is achieved. 

Closer to the right border of the dipolar window the 
front becomes turbulent and looking like precipitation, as 
seen in the lower image in Fig. [4] for B^ = 0.56 T. Given 
its sonic speed that cannot be provided by any thermal 
mechanism, one can speak of quantum detonation. The 
relaxation event in the figure began spontaneously on a 
surface near the left end where the resonance condition 
was satisfied; From there a regular slow-burning front is 
moving to the left and a quantum detonation front is 
going to the right. 

In thicker crystals, such diS = Ly = 1 mm in Fig. |5j 
the front speed is higher and non-flatness of the front is 
stronger pronounced. 

Although in the quantum detonation front the dipolar 
bias is wildly changing in space so that the system does 
not stick to the resonance (that was the main argument in 
the Id theory of fronts of tunneling [HI [13]), tunneling 
is still very strong to ensure sonic front speeds. This 
requires a closer investigation. Still, one can argue that 
dynamic crossing the resonance can facilitate transitions 
via Landau-Zener effect. 

To summarize, a full three-dimensional theory of 
dipolar-controlled quantum deflagration with the tem- 
perature effect, going over into quantum detonation, has 
been proposed. The quantum front can reach sonic 
speeds near biased ground-state tunneling resonances 
and is non-flat and turbulent. To activate quantum tran- 
sitions out of the metastable ground state, a strong trans- 



verse field has to be applied. An experimental challenge 
is to prepare the initial state that has a limited life time. 
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presentation. 
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